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Cold atom systems offer a great potential for the future design of new mesoscopic quantum systems with prop- 
erties that are fundamentally different from semiconductor nanostructures, such as quantum dots and quantum 
wires with electrons. Here, we investigate the analog of a quantum wire using ultra-cold particles, and find a 
new scenario for the quantum transport: Attractive interactions may lead to a complete suppression of current in 
the low-bias range, a total current blockade. We demonstrate this effect for the example of ultra-cold quantum 
gases with dipolar interactions. 



The electronic Coulomb blockade in mesoscopic quantum 
dots has been an intensive topic of research over the last 
two decades. The flow of electrons through a quantum dot 
between two reservoirs turned out to be an extremely versat- 
ile tool for addressing a wide range of fundamental effects. 
Examples range from investigating the structure of electronic 
many -particle states [1,2] and Kondo physics [3-5], to quan- 
tifying the spin dephasing due to coupling to nuclear degrees 
of freedom [6-8], or coherent effects [9]. 

Ultra-cold atoms in traps are very similar to quantum dots - 
a few quantum particles confined by a (often low-dimensional 
and harmonic) potential. What makes these systems partic- 
ularly interesting, is that one essentially can freely engineer 
their properties, and even control the shape and strength of 
the inter-particle interactions. This has sparked great interest 
in making devices with systems of ultra-cold atoms and mo- 
lecules analogous to those studied in electronics and spintron- 
ics [10-13]. 

In quantum-optical systems, initially the quantum confine- 
ment was achieved for very large systems with millions of 
atoms. Only recently, the few-body limit has been realized in 
a remarkable experiment by Serwane et al. [14]. They repor- 
ted the trapping of up to ten ^Li atoms where, intriguingly, 
they could determine the atom numbers down to single-atom 
precision, with full control over confinement and inter-particle 
interactions. 

"Interaction blockade" as the cold-atom analog of elec- 
tronic Coulomb blockade [15] was experimentally first seen in 
tunneling processes in optical lattices [16] and analyzed the- 
oretically for one-dimensional triple-well systems [17]. How- 
ever, the experimental realization of quantum transport of 
cold atoms through a small quantum few-body system that 
is brought in contact with two large atomic reservoirs has up 
to now posed a great experimental challenge. A very recent 
breakthrough by Brantut et al. [18], however, demonstrates the 
possibility to engineer both a ballistic and a diffusive channel 
between two cold atom reservoirs, opening up a host of new 
perspectives in mesoscopic quantum physics. 

Inspired by this recent experimental progress, here we study 
the quantum transport through quantum-dot-like confinement 
of a few ultra-cold fermions with an electric dipole moment. 
The possibility to align the dipole moment by an external field, 
allows tuning of the interaction [19]. This is shown to facilit- 
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Figure 1: Upper panel: Schematic figure of the system. Lower 
panel: Sketch of the setup in analogy to the case of mesoscopic con- 
ductors. Two reservoirs with a degenerate gas of ultra-cold spin-half 
dipolar particles are connected via a quasi one-dimensional structure, 
a "wire" of length 2a. The difference in chemical potential between 
the reservoirs, A/x, creates a particle current if the dipoles can be 
added and removed from the wire. Levels in the wire may be tuned 
by a gate potential, /Xgate- The interaction between the particles in the 
wire can be varied by the tilt angle of the dipoles, G, and allows to 
observe significantly different current patterns. 



ate studying localization effects and even offers a total current 
blockade, where the attractive interaction hinders transport for 
finite biases independent of the gate potential. 

The system considered in this Letter is sketched in Fig. 1. 
We assume the presence of two degenerate fermionic atomic 
clouds where the difference in chemical potential A/i can be 
controlled similar to the recent study in [18]. These clouds act 
as particle reservoirs, which are connected via a quasi one- 
dimensional confinement established by an appropriate trap. 
In this region, the potential energy of the particles can be var- 
ied by the parameter /igate in full analogy to electrons in gated 
semiconductor nanostructures. 

We consider a gas of particles with electric dipole moment, 
which can be orientated along an external field by tilt angle 
6 with respect to the z axis along the quasi one-dimensional 
channel (see Fig. 1). Within this confinement the interaction 
between the dipoles is of the form [20, 21] 

V!i{zr,Z2) = /7dd(6)Udd(|^i - Z2\/l±) + aS{z, - Z2) . (1) 
where lj_ is the characteristic length of the tight harmonic con- 
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finement in x and y. If the dipoles are aligned in the z direc- 
tion (9 = 0°) they attract each other, [/dd < 0, while they 
repel one another, Udd > 0, if they are orientated perpendic- 
ular to the z direction (6 = 90°). For an intermediate angle 
(Ocrit — 54.7°) this long-range part of the dipole interaction 
vanishes, i.e. Udd = 0. In addition to the common interaction 
between separated dipoles there is a contact term aS{zi — Z2) 
where a > for electric dipoles, as outlined in [21], where 
further details are given. 

For an arbitrary number of dipoles in the quasi one- 
dimensional confinement, we evaluate the quantum mechan- 
ical eigenstates using exact diagonalization. The transition 
between states of different particle numbers are due to trans- 
itions to the contacts with rates F [21]. This provides a Pauli 
master equation for the probabilities to find any many-particle 
state in the confined central region, which is solved numeric- 
ally similar to the case of electronic systems [22]. Thus we 
obtain the (particle) current N between the reservoirs and the 
(differential) conductance G = d7V/dA/i. 

Main results. — First we neglect the contact term of the di- 
polar interaction, assuming it to be eliminated by Feshbach 
resonances [23], and obtain the conductance diagrams dis- 
played in Fig. 2. For repulsive interaction between the dipoles, 
see Figs. 2(a)-(b), they resemble the Coulomb diamonds as 
intensively studied by electron transport in mesoscopic struc- 
tures [1]. This demonstrates the universality of the concept of 
interaction blockade. In these diamonds of transport block- 
ade (areas labeled by TV = 1, 2, 3, 4, 5) the current between 
the reservoirs is strongly suppressed. In contrast to electronic 
systems, the tunability of the interaction for dipolar fermions 
allows to reduce the interaction strength. Figs. 2(c)-(d), and 
even reach a scenario where the interactions become attract- 
ive, see Fig. 2(e). 

In the last case. Fig. 2(e), we obtain a total current blockade 
at low detuning A/i independent of the gate potential /igate- 
This is a clear-cut signature of attractive interaction. 

In general, these features can be understood by the two- 
fold degeneracy of the single particle levels due to the particle 
spin. For the conventional repulsive Coulomb interaction, the 
situation of a single particle in a shell is stable, as adding a 
second particle requires the charging energy U > 0. Thus, 
one observes the Coulomb diamonds with an odd number of 
particles N and lines of finite conductance at the separation 
to the Coulomb diamonds with even N (Figs. 2(a), (b)). With 
decreasing interaction Udd, the width of all diamonds shrinks 
and the width of the odd- A/" diamonds vanishes at /7dd = as 
can be seen in Fig. 2(d). Now, for negative Udd the situation 
of a single fermion in a shell is unstable as it immediately at- 
tracts a particle with the opposite spin. This instability does 
not allow for configurations with odd N for low A/i. There- 
fore, single-particle transitions between the reservoir and the 
central region are excluded, resulting in the absence of current 
flow in the region of total current blockade, see the magenta 
shaded area in Fig. 2(e). (The case of two-particle transitions 
is addressed below.) 

For weak interactions the structure of the conductance 
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Figure 2: Conductance between the particle reservoirs as a function 
of reservoir potential difference A/x and gate potential /Xgate- Here the 
contact part of the dipolar interaction is neglected and the long-range 
part, which can be tuned by the angle G of the external field, changes 
from (a) strong repulsive, via (d) non-interacting, to (e) the weak 
attractive case. The region of total current blockade for attractive 
interaction is colored in magenta in (e). The red lines in (c) and (e) 
indicate the results of a simplified quasi-independent-particle model. 
The /i-scale is in units of Ti^ /mo? . Note the different scales in panel 
(a) and (e). 

peaks in Figs. 2(c), (e) can be understood within a standard 
quasi-independent-particle model, as outlined in [21]. The 
corresponding red lines provide results similar to the main 
conductance lines obtained from the full many -particle calcu- 
lation. Thus correlations do not play any essential role here. 

In contrast, such an approach does not hold for stronger in- 
teraction strengths. Here the many -particle states show strong 
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Figure 3: Particle density (left) and pair-correlation function (right) 
for = 2 particles at the tilt angles G used in Fig. 2(a)-(e). For the 
pair-correlation function one particle is fixed at the position marked 
with the symbol 0. As the interaction goes from strong repulsive 
(G = 90°) to weakly attractive (G = 54.2°) the two particles evolve 
from a localized state to a delocalized state with a slight tendency to 
clustering. 



localization effects as shown in Fig. 3 for the two-particle 
states. For 6 = 90° and to a smaller extent for 6 = 60°, one 
observes two peaks in the particle density (left panel), and the 
pair-correlation function (right panel) shows that the probab- 
ility to find the two particles within the same peak is strongly 
reduced. This is the scenario of Wigner localization as very 
recently studied theoretically for cold polar molecules in [24] . 
In full analogy to mesoscopic electron conduction [22], sig- 
natures of this localization can be clearly detected in the con- 
ductance plots Fig. 2(a)-(b) where several, almost degenerate 
lines are observed on the top of the diamonds, which result 
from spin excitations of the localized particles. 

Improved interaction model. — For attractive interaction 
(9 = 54.2°), the pair-correlation function is shifted to the 
right, see the right panel of Fig. 3, i.e. the probability to find 
both fermions on the same spot is enhanced for the ground- 
state. In this case the contact interaction in Eq. (1) is of 
particular relevance. Taking this term into account provides 
some modifications of the scenario depicted in Fig. 2, while 
the main features remain. Here we study the case of electric 
dipoles, where the contact interaction is repulsive {ol > 0). 
The additional repulsion compensates a part of the long-range 
attraction, so that smaller angles 9 are required to observe the 
vanishing of the diamonds with odd N. Furtermore, since the 
particle density increases with the number of particles TV, the 
contact interaction becomes more relevant for higher N, and 
thus smaller angles are required for the vanishing of diamonds 
with higher N. We have observed this for e.g. 9 = 46°, 
where the A/" = 1 diamond has already vanished, while the 
= 3 diamond is very small and the = 5 diamond is 
still well established. In this case the total current blockade 
due to the attractive interaction extends only over a part of the 
spectrum. 

Pair-tunneling. — As discussed above, the situation of a 
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Figure 4: Conductance through a single spin-degenerate level at 
energy Sd for the case of inter-particle attraction, resulting in negative 
charging energy U < 0. The temperature is ksT = |[/|/10, while 
the couplings are given by = Fr = |[/|/50, and Fl = F^ = 
I [/ 1/4 for (a) and (b), respectively. The red and blue lines show the 
onset of sequential and pair- tunneling, respectively. 



single fermion in a shell is unstable for the case of attracting 
particles, U^d < 0. Thus, single-particle transitions between 
the reservoir and the central region are excluded for suffi- 
ciently low vales of temperature and bias A/i. Here we want 
to illuminate the role of two-particle transitions, which may 
occur due to higher-order processes in the coupling between 
the reservoirs and the central region [25, 26]. 

There are two kinds of processes: Normal co-tunneling, 
which results in a weak background conductance for any bias, 
and pair-tunneling which, neglecting the effects of temperat- 
ure and lifetime broadening, only is allowed for \E2n-\-2 — 
E2n\ < A/i/ 2, where En '^^ the ground state energy of 
the A^-particle state and n is an integer. When present, 
pair-tunneling gives a more pronounced contribution than co- 
tunneling [25]. Being of second order, these processes scale 
as F^, where F is the single particle transition rate. Thus, for 
weak enough couplings, they can be neglected compared to 
sequential single-particle tunneling. 

Figure 4 shows the differential conductance of a single spin- 
degenerate level with negative charging energy /7 < 0, calcu- 
lated by the second order von Neumann formalism [27, 28]. 
For weak contact coupling. Fig. 4(a) displays only a small 
conductance at low values of A/i. This can be attributed 
to a weak pair-tunneling background and to the temperature 
broadening ~ of the direct tunneling peaks at Sd = 

-U/2 ± {Ajj. + U) /2 for A/i > -[/, which correspond to the 
red lines in Fig. 4. This demonstrates that the total blockade 
of conductance is verified for F <C /c^T <C |/7|, as is the case 
in Fig. 2(e). 

On the other hand, one has to keep in mind that, as F ap- 
proaches U, pair-tunneling becomes energetically allowed. 
Hence, we observe the onset of conduction along the blue 
lines \E2—Eo\ = A/i in Fig. 4(b). (In our case E2 = 2s d + U 
and £^0 = 0.) Normal co-tunneling can also be observed as 
a weak backgrond present at all A/i and £d- Thus, the total 
blockade of conductance does not persist at strong couplings 
between the central region and the reservoirs. For even higher 
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couplings, our model fails and Kondo-like effects become im- 
portant [29]. This shows that a certain confinement of the 
central region is required for the observation of the total cur- 
rent blockade as otherwise pair-tunneling masks the scenario. 

Experimental challenges. — From the experimental point of 
view, measuring a weak atomic current in a mesoscopic trans- 
port process appears challenging. First experimental studies 
of quantum transport through atom traps have very recently 
been reported by Brantut et al. [18], where integrated current 
is measured by a sensitive detection of population differences 
in the reservoirs. This opens up a new field of mesoscopic 
physics research. Complementary experimental information 
on the atomic current could, for instance, be inferred from 
a time-of-flight absorption image that renders the momentum 
distribution of the transported atoms. As an alternative, a stim- 
ulated Raman adiabatic passage (STIRAP) of the atoms could 
be induced by irradiating the transport region with two spa- 
tially displaced laser beams (see, e.g., Ref. [30]). An atom that 
propagates through this irradiated region would then necessar- 
ily transfer a photon from one of the laser beams to the other, 
while an atom that propagates in the opposite direction would 
revert this photon transfer. A careful measurement of the net 
photon transfer between the beams after a suitable evolution 
time would then give rise to the integrated atomic net current 
across the atom-photon interaction region. We remark that 
standard techniques to detect individual atoms using fluores- 
cence imaging [31, 32] or electron beams [33] would not work 
in this context as they do not distinguish between left-moving 
and right-moving atoms. 

Conclusions. — We have shown that dipolar quantum gases 
allow for the observation of a total current blockade for small 
differences in chemical potentials between the reservoirs. In 
this context the often neglected contact interaction part of the 
dipole-dipole interaction turns out to repress the onset of total 
current blockade. 

From the experimental side, studies of quantum transport 
with ultra-cold atoms and the many-body effects of interaction 
blockade are still in their infancy. Here, we highlighted the 
prospects for the specific example of a few-body system with 
dipolar interactions between the confined atoms. We demon- 
strated the possibilities offered by the tunability of the dipole- 
dipole interaction in a quasi one-dimensional geometry by an 
external field. 
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THEORETICAL MODELLING 

The interaction between two ideal polarized dipoles may be 
expressed as 

14d = ^ (1 - 3 cos2 0,,) + ^Cd^S\r) (2) 

where r is the distance between the dipoles, and O^d is the 
angle between the vector r connecting the dipoles and the 
vector d pointing along the dipole moment [2]. For electric 
dipoles the coupling strength is cP = p^/(47reo), where p is 
the dipole moment strength, eo is the vacuum permittivity and 
C = 1, while for magnetic dipoles the coupling strength is 
d'^ = iiQg\ii\/ {^it), where /io is the vacuum permeability, 
ql the Lande factor, (Ib the Bohr magneton and C = —2. 
While the first term provides the common angular depend- 
ence of dipole-dipole interaction, the second term provides a 
contact interaction which is frequently disregarded. 

The dipoles are confined in x and ^ by a 2-dimensional har- 
monic oscillator of characteristic length rendering a quasi 
one-dimensional system in the z-direction for small l^. Integ- 
rating over the lateral x and y degrees of freedom one arrives 
at an effective one-dimensional dipole-dipole interaction 

— 2 Cd'^ 

V!!{z) = UM\z\/l^) + ^^5{z) . (3) 

Here 

d^[l + 3cos(29)] 
= 8Zi ' 

where 6 is the dipole tilt angle, i.e. the angle between the 
vector d and the z-axis, and 

VddH = -2u + V^(l + ii^)e^'/^ erfc(ii/V2) (5) 

where erfc is the complementary error function, as shown 
in [3]. Note that depending on the dipole tilt angle B, the 
interaction coefficient /7dd can be either positive or negat- 
ive. Hence, by changing 6, the interaction can be tuned to 
be repulsive, zero or negative, its strength /7dd ranging from 
-d^/2l\ to d^/Al\. 

In the z-direction the wire is modelled as a finite square 
well (see Fig. 1 in the main text) of width 2a and barrier 
height ^0- Applying the single-particle basis of eigenstates 
for this potential well, the configuration interaction method 
(exact diagonalization) is used to find the lowest energy states 



of = 1 to = 6 dipoles in the quantum wire. Here 
the dipolar particles are assumed to be spin-half fermions. 
The parameters used were d^ = Ih^a/m, l± = 0.14a and 
Vo = 300h^/ma'^. 

Transition coefficients Ti)a{kcri) from the wire to state k 
in the contacts are found from creation overlaps {b\d^a) 
between many particle states \a) and \ b) with spin a, evaluated 
in the contact to lead i following the work of [4, 5] for meso- 
scopic electric systems. The transition rates are then given by 
Fermi's golden rule 

Ta^b/ = yY1 \Tba{kcTi)\^S{E, -Ea- E^) . (6) 

k(j 

Neglecting the off-diagonal terms of the density matrix a 
Pauli master equation can be written in terms of these trans- 
ition rates. In our calculations, the temperature was set to 
ksT = {)mn^/ma^. 

A simple analytical expression for the ground state levels 
can be obtained within an quasi-independent-particle shell 
model. The one-particle energies of the quantum well are 
approximated by n^Ei where n = 1,2,... and Ei is the 
one-paricle ground state energy, obtained numerically. Using 
the analytic eigenf unction of the infinite well, we approximate 
the interaction energy by first order perturbation theory. Then 
the energy difference between the A" + 1 and the A^-particle 
ground state is simply given by 

/iAT+l = /igate + + if E^ + (n + 5)Uq (7) 

where Uq = 2l±Udd/ci and /igate is the gate potential relative 
to the bottom of the well. Here n = [N/2\ and = for 
even N = 2n and ^ = 1.5 for odd A" = 2n + 1. The lines of 
the diamonds are given by the crossing points of /iAr+i with 
the chemical potential ±A/i/2 in the left or right reservoir, 
respectively, resulting in condition //at+i = ±A/i/2 for the 
dashed red lines displayed in Fig. 2 in the main text. 
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